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Abstract 

Background: Recent studies have clearly demonstrated the enormous virus diversity that exists among wild 
animals. This exemplifies the required expansion of our knowledge of the virus diversity present in wildlife, as well 
as the potential transmission of these viruses to domestic animals or humans. 

Methods: In the present study we evaluated the viral diversity of fecal samples (n = 42) collected from 10 different 
species of wild small carnivores inhabiting the northern part of Spain using random PCR in combination with 
next-generation sequencing. Samples were collected from American mink (Neovison vlson), European mink 
{Mustela lutreola), European polecat {Mustela putorius), European pine marten {Martes martes), stone marten 
(Martes foina), Eurasian otter {Lutra lutra) and Eurasian badger {Meles meles) of the family of Mustelidae; 
common genet (Genetta genetta) of the family of Viverridae; red fox {Vulpes vulpes) of the family of Canidae 
and European wild cat {Felis silvestris) of the family of Felidae. 

Results: A number of sequences of possible novel viruses or virus variants were detected, including a theilovirus, 
phleboviruses, an amdovirus, a kobuvirus and picobirnaviruses. 

Conclusions: Using random PCR in combination with next generation sequencing, sequences of various novel viruses 
or virus variants were detected in fecal samples collected from Spanish carnivores. Detected novel viruses highlight the 
viral diversity that is present in fecal material of wild carnivores. 

Keywords: Theilovirus, Viral metagenomics, Picobirnavirus, Amdovirus, Phlebovirus, Carnivore, Wild carnivore, Wildlife, 
Red fox, Fox, Otter, Common gennet, Bunyavirus 



Background 

Transmission of viruses from wildlife to humans con- 
tinues to cause outbreaks of disease in humans. Exam- 
ples of recent outbreaks are the Middle East Respiratory 
Syndrome-coronavirus (MERS-CoV) that may have orig- 
inated from bats and/or camelids and the influenza A 
(H7N9) virus that originated from wild birds [1-5]. A 
systematic exploration of viruses present in several key 
host species of wild animals might provide important in- 
formation to find the original host or carriers of viruses 
of future outbreaks of viral disease among domestic ani- 
mals, endangered animal species, and humans [6]. Fur- 
thermore, information about the presence of viruses in 
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healthy hosts provides a baseline level for viruses present 
in these animals in case an outbreak of disease occurs. 
In previous viral metagenomics studies, high numbers of 
new viruses have been identified [7-10]. The results of 
these studies have highlighted that our knowledge of the 
viral reservoir is far from complete and many, as yet, un- 
identified viruses circulate among humans and wild and 
domestic animals. However, there is an enormous diver- 
sity of viral sequences and viral metagenomics efforts 
should be focused on outbreaks of disease and viral meta- 
genomics on samples collected from a selected number of 
key species [11]. 

Wild carnivores are known carriers of several viral 
pathogens that can affect domestic animals and humans, 
including rabies and canine distemper virus [12,13]. In 
addition, in previous studies various previously unknown 
viruses have been detected in European badgers, red 
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foxes and European pine martens in the Netherlands 
[9,14]. In the present study, we evaluated the viral di- 
versity of fecal swabs or fecal specimens collected from 
10 different small carnivore species of the Mustelidae, 
Canidae, Viverridae and Felidae families inhabiting 
northern Spain. 

Results and discussion 

Metagenomic overview 

Using random amplification in combination with next- 
generation sequencing, more than 320,000 trimmed se- 
quence reads were obtained of fecal samples collected 
from the carnivores of the present study (Figure 1). Reads 
were classified into eukaryotic viruses, phages, bacteria 



and eukaryotes. Many of the identified sequences were of 
bacterial or eukaryotic origin. A substantial proportion of 
the reads did not have any significant hits for nucleotide 
or amino acid sequences in GenBank. In addition, several 
reads were detected that had the closest similarity to vi- 
ruses. In the majority of the samples, sequences of the 
order Caudovimles were detected and in 26 out of 42 
samples, sequences were detected that had the closest 
similarity to viruses known to infect eukaryotes (Figure 2A, 
Table 1). Viruses belonging to the families of Anelloviridae, 
Astroviridae, Bunyaviridae, Caliciviridae, Circoviridiae, 
Parvoviridae subfamily Parvovirinae, Picobirnaviridae, 
Picornaviridae, Rhabdoviridae, and Retroviridae were 
detected (Figure 2B). Furthermore, sequences were detected 




ANALYSED SPECIMENS 
Mustelidae 

■ European mnk [Mustela kAreoia) n = 9 
® Amencan mink {Neovison vtson) n = 10 

• European polecat [Musteia putonusi n ■ 3 



Canidae 

▲ Red fox (Vuipes wipes) n-4 

Fel.flae 

V European wtfdcal {Felrs sitvestns) n - 1 
Vivc"'dacc 

% Common genet ( Genefla genetta) n ■ 7 



Highways and National roads 
Uroan Areas 



0 European pine marten (Maries manes) n = 2 
0 Stone marten {Maries fotna) n = 1 

0 Eurasian badger (Metes metes) n = 4 

0 Eurasian otter {Lutra lutra) n = 1 

Figure 1 Map with locations of sample collection. Maps showing the location of the study area in the Basque Country and La Rioja regions 
(Spain) and spatial distribution of the 42 analysed carnivore samples of the Mustelidae, Canidae, Felidae and Viverridae families. 
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American mink - 
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Anellovirus 
Circo-like virus 
Circovirus 
Cyclovirus 
Meles meles fecal virus 
NIH-CQV 
Papillomavirus 
Amdovirus 
Dependovirus 
Parvovirus - 
Astrovirus 
Phlebovirus 
Vesivirus 
Rhabdovirus 
Picobirnavirus 
Tetnovirus 
Noravirus 
Cardiovirus 
Kobuvirus- 
Rosavirus 
Betaretrovirus 
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Number of samples with viral sequences 



lAmencan mink 
]Common genet 
I Eurasian badger 
I Eurasian otter 
JEuropean mink 



I European pine marten 
I European polecat 
]Red fox 
I Stone marten 
JWild cat 



Figure 2 Overview of detected viral sequences in the present 
study. Number of animals included in this study (bars) and the 
number of samples in which viruses where detected (grey areas; 
A). Number of animals in which viral sequences were detected for 
each genus/ group (B). 



that had the closest similarity to the recently proposed fam- 
ily of Breviviridae and the recently described hybrid DNA 
virus NIH-CQV/PHV which was identified as a contamin- 
ant of silica column-based nucleic acid extraction kits 
[9,15,16]. No sequences were detected that were identical 
to currently known zoonotic viruses. A proportion of the 
detected viral sequences had the closest similarity to viruses 
previously detected in birds and rodents. For example, in 
an European mink (sample 26), sequences were detected 



with >95% homology on the nucleotide level with Turkey 
parvovirus and in a stone marten (sample 41), sequences 
were detected with 94-96% homology on the nucleotide 
level with Encephalomyocarditis virus type 2 isolate RD 
1338 (D28/05) detected in a wood mouse (Apodemus sylva- 
ticus) [17]. These viruses most likely originate from the diet 
of the animals. In addition, sequences with >95% identity 
on the nucleotide level to viruses that are known to infect 
mink were detected in European and American mink, in- 
cluding Mink calicivirus strain MCV-DL/2007/CN [18] 
(samples 1 and 8) and Aleutian mink disease virus (sample 
30). Antibodies to Aleutian mink disease parvovirus have 
been detected in a cohort of free-ranging European mink in 
southwestern France and northern Spain previously, but 
not in another cohort of free-ranging European mink in 
Navarra, Spain [19-21]. Additional sampling and confirm- 
ation by specific PCR is necessary to indeed confirm that 
the Aleutian mink disease parvovirus is circulating among 
these animals. Besides these sequences that had high hom- 
ology with known viruses, also sequences were detected 
that had the closest similarity to viruses, but with only low 
homology. A number of sequences of potentially novel vi- 
ruses or virus variants, including a theilovirus, phlebo- 
viruses, an amdovirus, a kobuvirus and picobirnaviruses, 
were further characterized in the present manuscript, while 
sequences of the other viruses are preliminary and need 
further characterization. 

Genet fecal theilovirus 

Picornaviruses are small, positive-stranded non-enveloped 
RNA viruses. This large family was subdivided into several 
genera [22]. The genus Cardiovirus of the family Picorna- 
viridae currently consists of two species, Theilovirus and 
Encephalomyocarditis virus (EMCV) [22] . Viruses belong- 
ing to the species Theilovirus were detected initially 
mainly in rodents. In these animals, viruses cause primar- 
ily infection of the digestive tract without clinical signs, 
but extra-intestinal infection occurs and can cause an 
acute encephalomyelitis and a chronic demyelinating in- 
fection of the central nervous system [23]. Vilyuisk hu- 
man encephalomyelitis virus, another strain of the 
species Theilovirus, was isolated from humans with en- 
cephalomyelitis, but only after serial passage over mice 
brains [24,25]. More recently, Saffold virus was discov- 
ered in a stool sample of a child with fever of unknown 
origin, and additional research revealed that infection 
with this virus was common [26-28] . 

We detected in fecal material of a stone marten 
(sample 41) sequences with high similarity to a novel 
Encephalomyocarditis type 2 virus detected recently in 
a wood mouse {Apodemus sylvaticus) [17]. In addition, 
sequences were detected that had the closest similarity 
to viruses of the species Theilovirus in a rectal swab col- 
lected from a common genet (sample 14). The obtained 



Bodewes et al. Virology Journal 2014, 11:89 
http://www.virologyj.eom/content/1 1/1/89 



Page 4 of 13 



Table 1 Overview of samples used in this study and detected mammalian viral sequences 



No 


Species 


Date of collection 

^uu/ 1 1 1 1 1 1/ y y y y / 


Observations 1 


Analyzed 2 


Reads 


Detected viral 

3 4 

JCUUCI ILej 


1 


American mink [Neovison vison) 


08/1 1/2012 


e 


F 


5316 


t.b.c. 


2 


American mink [Neovison vison) 


15/1 1/2012 


e 


RS 


4168 


t.b.c. 


3 


American mink [Neovison vison) 


10/1 1/2012 


e 


F 


7563 


t.b.c. 


4 


American mink (Neovison vison) 


09/1 1/2012 


e 


F 


9041 


n.d. 


5 


American mink (Neovison vison) 


19/09/2012 


e 


F 


4919 


t.b.c. 


6 


American mink [Neovison vison) 


13/10/2012 


e 


F 


4454 


n.d. 


7 


American mink [Neovison vison) 


21/09/2012 


e 


F 


4587 


n.d. 


8 


American mink [Neovison vison) 


26/10/2012 


e 


F 


5592 


t.b.c. 


9 


American mink [Neovison vison) 


1 1/09/2012 


e 


F 


5838 


n.d. 


10 


American mink [Neovison vison) 


10/04/2012 


e 


F 


6925 


n.d. 


1 1 


Common genet [Genetto genetto) 


08/1 1/2012 


It 


RS 


6353 


t.b.c. 


12 


Common genet [Genetto genetto) 


05/1 1/2012 


fd 


F 


8642 


t.b.c. 


13 


Common genet [Genetto genetto) 


04/07/2012 


It 


RS 


5582 


t.b.c. 


14 


Common genet [Genetto genetto) 


09/1 1/2012 


fd 


F 


5217 


Cardiovirus, Picobirnavirus, t.b.c. 


15 


Common genet [Genetto genetto) 


14/12/2012 


fd 


RS 


6758 


n.d. 


16 


Common genet [Genetto genetto) 


19/02/2013 


fd 


RS 


7327 


t.b.c. 


17 


Common genet [Genetto genetto) 


10/05/2013 


fd 


RS 


9068 


t.b.c. 


18 


Eurasian badger [Meles Meies) 


04/01/2013 


fd 


RS 


6465 


t.b.c. 


19 


Eurasian badger [Meles Meies) 


13/05/2013 


fd 


RS 


1 1228 


n.d. 


20 


Eurasian badger [Meles Meles) 


12/05/2013 


fd 


RS 


3763 


n.d. 


21 


Eurasian badger [Meles Meles) 


21/05/2013 


fd 


RS 


6946 


n.d. 


22 


Eurasian otter (Lutro lutro) 


09/01/2013 


fd 


F 


5500 


Phlebovirus, t.b.c. 


23 


European mink [Mustelo lutreolo) 


06/1 1/2012 


It 


RS 


1 1356 


t.b.c. 


24 


European mink [Mustelo lutreolo) 


07/12/2012 


It 


RS 


13933 


t.b.c. 


25 


European mink [Mustelo lutreolo) 


08/1 1/2012 


It 


RS 


13119 


n.d. 


26 


European mink [Mustelo lutreolo) 


16/1 1/2012 


It 


RS 


7242 


t.b.c. 


27 


European mink [Mustelo lutreolo) 


26/1 1/2012 


fd 


RS 


10046 


t.b.c. 


28 


European mink [Mustelo lutreolo) 


7/03/2013 


It 


RS 


9027 


n.d. 


29 


European mink [Mustelo lutreolo) 


8/03/2013 


It 


RS 


7391 


n.d. 


30 


European mink [Mustelo lutreolo) 


6/03/2013 


It 


RS 


6593 


t.b.c. 


31 


European mink [Mustelo lutreolo) 


20/03/2012 


fd 


RS 


9418 


n.d. 


32 


European pine marten [Mortes mortes) 


11/1 1/2012 


fd 


F 


9624 


n.d. 


33 


European pine marten [Mortes mortes) 


05/03/2013 


fd 


RS 


9961 


t.b.c. 


34 


European polecat [Mustelo putorius) 


19/02/2013 


It 


RS 


7422 


t.b.c. 


35 


European polecat [Mustelo putorius) 


zo/Uz/zU 1 J 


fA 


DC 


9741 


n.a. 


36 


European polecat [Mustelo putorius) 


15/05/2013 


fd 


RS 


17970 


n.d. 


37 


Red fox [Vulpes vulpes) 


08/01/2013 


fd 


RS 


5163 


Kobuvirus, Phlebovirus 


38 


Red fox [Vulpes vulpes) 


29/01/2013 


fd 


RS 


5722 


n.d. 


39 


Red fox [Vulpes vulpes) 


17/03/2013 


fd 


F 


5842 


Kobuvirus, t.b.c. 


40 


Red fox [Vulpes vulpes) 


08/04/2013 


fd 


F 


8828 


Picobirnavirus, Amdovirus, t.b.c. 


41 


Stone marten (Martes foina) 


15/11/2012 


fd 


F 


6809 


t.b.c. 


42 


Wild cat [Felis sllvestris) 


24/05/2013 


fd 


RS 


7224 


t.b.c. 



1 e: euthanized animal (in the framework of an eradication program of invasive species), It: live trapped animal, fd: found dead. 
2 F: fecal sample, R: rectal swab. 
3 nd: no viral sequences detected. 

4 t.b.c: one or more sequences were detected that were most closely related to viral sequences, but sequences were not characterized in the present manuscript. 
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sequences were further analyzed by designing primers 
based on the obtained sequence reads and the partial 
sequence (6488 nt) of a novel theilovirus, tentatively 
called Genet fecal theilovirus (GFTV; Genbank acces- 
sion KF823815), was obtained. Despite various attempts 
using degenerate primers and 3' RACE PCRs, no more 
sequences could be obtained of the 5' and 3' end, pos- 
sibly due to low virus genome copy numbers present in 
the original material. The obtained sequence contained 
the partial non-coding region of the 5'end and the par- 
tial ORF encoding for the polyprotein gene. Based on 
splice-site predictions and alignment with other theilo- 
viruses, the complete coding sequence of the L protein 
(228 nt), the complete PI (2462 nt) capsid-encoding re- 
gions, the complete P2 (1785 nt) non-structural coding 
region and the partial P3 (1455 nt) non-structural cod- 
ing region were detected in the partial polyprotein. Pair- 
wise identity analysis and phylogenetic analysis of the 
nucleotide sequences of the complete Leader gene, PI, 
P2 and partial P3 gene showed that Genet fecal theilovirus 
probably belongs to a novel genotype, with maximum 
pairwise identities on the nucleotide (and deduced amino 
acid) level of respectively 69 (70), 63 (68), 62 (65) and 82 
(65)% (Figure 3). Additional alignment of the deduced 
amino acid sequence of the PI gene with the major sur- 
face structures of other viruses of the species theilovirus 
(VP2 puffs A and B, the VP3 knob and VP1 loops 1 and 2) 
indeed shows the high divergence between the potential 
immunogenic sites of these viruses (Additional file 1: 
Figure SI) [29,30]. In addition to the major ORF encoding 
the polyprotein, an alternative ORF of 423 nt (140AA) 
was present in the Genet fecal theilovirus sequence based 
on an alternative initiation codon 13 nt downstream of 
the authentic initiation codon, which indicates that the L* 
protein previously observed in Theiler's murine encepha- 
lomyocarditis virus and other strains of the species Theilo- 
virus is also present in this novel virus [29,31]. Since this 
animal was found dead, possibly due to a car accident, 
and the carcass was stored at -20 for a few months, a nec- 
ropsy could be performed. The carcass was defrosted but 
no abnormalities were detected in this animal upon 
macroscopic and microscopic examination of various tis- 
sues, including the brain. This indicates that the novel 
Genet fecal theilovirus has not caused any significant dis- 
ease. Since picornaviruses are known to be very resistant 
to a low pH, further research needs to be performed to 
elucidate whether the detected virus is derived from the 
prey species of this animal or is a novel theilovirus of the 
common genet. 

Novel phleboviruses in feces of a red fox and an eurasian 
otter 

Bunyaviruses are segmented, negative-sense single stranded 
RNA viruses. Genomes consist of three different segments, 



the Large (L), Medium (M) and Small (S) segment. At 
present the ICTV has recognized five different genera, 
Orthobunyavirus, Hantavirus, Nairovirus, Phlebovirus, 
and Tospovirus [22]. The genus Phlebovirus consists of a 
genetically diverse group of viruses, some of which were 
described very recently [32-34]. A few members of this 
genus were identified as important pathogens in humans 
and domestic animals, including Rift valley fever virus and 
the Severe fever with thrombocytopenia syndrome virus 
or Huaiyangshan virus [35,36]. 

In fecal material of an Eurasian otter (sample 22) and 
a red fox (sample 37) sequences were detected that had 
the closest similarity to viruses of the genus Phlebovirus 
(Table 1). Additional virus reads were obtained of the 
samples of these animals (total 43,584 of sample 22 and 
30,064 reads of sample 37). By analysis of the additional 
obtained sequences of sample 22, the partial viral nucleo- 
protein (NP) gene (S segment) and partial glycoprotein 
(G) gene (M segment) of a novel phlebovirus, tentatively 
called Otter fecal phlebovirus (NP gene 699 nt, covered by 
16 reads, G gene 984 nt, covered by 59 reads, Genbank ac- 
cessions KF823816 and KF823817), were identified. In 
addition, by analysis of the additional obtained sequences 
of sample 37, the partial viral nucleoprotein (NP) gene 
(S segment) and partial glycoprotein (G) gene (M seg- 
ment) of another novel phlebovirus, tentatively called 
Red fox fecal phlebovirus (NP gene: 606 nt, covered by 
33 reads, G gene: 1110 nt, covered by 63 reads, Genbank 
accessions KF823818 and KF823819) were identified 
(Figure 4A, B). Also sequences were detected with the 
closest similarity to sequences of the large segment of 
phleboviruses, but with a lower coverage of reads. Pair- 
wise identity and phylogenetic analysis of the deduced 
amino acid sequence of the partial NP and Gl genes 
with various other viruses of the genus Phlebovirus sug- 
gest that both viruses are highly divergent viruses be- 
longing to the genus Phlebovirus with pairwise identities 
of less than 34% on the deduced amino acid level of 
both viruses of both genes (Additional file 1: Table SI, 
Table S2, Figure 4C, D). Of interest, phleboviruses are 
transmitted by arthropods or ticks, but these viruses 
were detected in the fecal material of two animals in this 
study. It has been demonstrated for Rift valley fever 
virus that after systemic infection virus can be detected 
in feces [37], therefore also these animals might have 
been infected systemically. However, the identification 
of these viruses in fecal content could be also due to the 
ingestion of preys (e.g. small mammals) with arthropods 
or ticks. 

Red fox fecal amdovirus 

Parvoviruses are small non-enveloped single-stranded 
DNA viruses. The family Parvoviridae has been subdi- 
vided into two different subfamilies, Parvovirinae and 



Bodewes et al. Virology Journal 2014, 11:89 
http://www.virologyj.eom/content/1 1/1/89 



Page 6 of 1 3 



Leader gene 




EU718732 TMEV 
EU723238TMEV 
M16020TMEV 
M20562 TMEV 
EU718733 TMEV 
M94868 VHEV 
M20301 TMEV 
HQ652539 TMEV 
AB090161 TRV 
EU542581 TRV 
%'EU815052TRV 
GFTV 
EU681 177 SAFV-2 
JF813004 SAFV-2 
GU943514 SAFV-3 
EF165067 SAFV-1 
HM181996 SAFV-3 
FJ463615SAFV-5 
FR682076 SAFV-2 
HM181999 SAFV-3 
GU943515 SAFV-3 
M81S61 EMCV-1 
AJ617362EMCV-1 
AY296731 EMCV-1 
M22457 EMCV-1 
122089 Mengo 
JX257003 EMCV-2 
AJ61 7358 EMCV-1 
DQ635185 EMCV-1 



P2 gene 




EU723238 TMEV 
EU718732 TMEV 
M 16020 TMEV 
M20301 TMEV 
HC6S2S39 TMEV 
U20288 TMEV 
EU718733 TMEV 
M20562 TMEV 
AB090161 TRV 
EU542S81 TRV 
EU815052TRV 
JF813004 SAFV-2 
FJ463615SAFV-5 
FR682076 SAFV-2 
EU681 1 77 SAFV-2 
HM181996 SAFV-3 
EF165067 SAFV-1 
HM181999 SAFV-3 
GU943514 SAFV-3 
•GFTV 

JX257003 EMCV-2 
M22457 EMCV-1 
AY296731 EMCV-1 
M81861 EMCV-1 
L22089 Mengo 




P1 gene 



EU718732 TMEV 
EU723238TMEV 
M 16020 TMEV 
M20562 TMEV 

H0652539 TMEV 
U20288 TMEV 
EU718733 TMEV 
M20301 TMEV 
-M94868 VHEV 
AB090161 TRV 
EU542581 TRV 
EU815052TRV 
HM181999 SAFV-3 
GU943515 SAFV-3 
GU943514 SAFV-3 
HM181 996 SAFV-3 
FJ463615SAFV-5 
EF165067 SAFV-1 
FR682076 SAFV-2 
EU681177 SAFV-2 
JF813004 SAFV-2 
GFTV 

JX2S7003 EMCV-2 
L22089 Mengo 
DQ83S1 85 EMCV-1 
AJ61 7358 EMCV-1 
M22457 EMCV-1 
AY296731 EMCV-1 
M81861 EMCV-1 
AJ617362 EMCV-1 



P3 gene 



DQ835185 EMCV-1 




EU723238 TMEV 
M20301 TMEV 
M20562 TMEV 
HCJ652539 TMEV 
EU718732 TMEV 

6020 TMEV 
EU718733 TMEV 
EU81S052TRV 
AB090161 TRV 
EU542581 TRV 
EU681 177 SAFV-2 
FR682076 SAFV-2 
FJ463615SAFV-5 
HM181996 SAFV-3 
EF165067 SAFV-1 
GU94351 4 SAFV-3 
JF813004 SAFV-2 
HM181999 SAFV-3 
-GFTV 

— JX257003 EMCV-2 

— L22089 Mengo 
-DQ8351 85 EMCV-1 
•M22457 EMCV-1 
-M81861 EMCV-1 
-AY296731 EMCV-1 



Figure 3 Phylogenetic analysis of Genet fecal theilovirus (GFTV). Phylogenetic neighbor-joining tree with p-distance and 1,000 bootstrap 
replicates of the nucleotide sequence of the Leader gene (A), PI gene (B), P2 gene (C) and partial P3 gene (D) of Genet fecal theiloviruses and 
various other viruses of the genus Cardiovirus. TMEV: Theiler's murine ecephalomyelitis vrus, VHEV: Vilyuisk human encephalomyelitis virus, TRV: 
Theiler's-like rat virus, SAFV: Saffold virus, EMCV: encephalomyocarditis virus. Indicated are genbank accession numbers and bootstrap values. 



Densovirinae. Viruses of the Densovirinae infect arthro- 
pods, while viruses of the Parvovirinae infect vertebrates. 
At present, the International Committee on Taxonomy of 
Viruses (ICTV) has recognized eight different genera of 
the subfamily Parvovirinae: Amdoparvovirus, Aveparvo- 
virus, Bocaparvovirus, Copiparvovirus, Dependoparvovirus, 



Erythroparvovirus, Protoparvovirus and Tetraparvovirus 
[22]. The genus Amdoparvovirus currently has two mem- 
bers, Aleutian mink disease virus and Gray fox amdovirus. 
Aleutian mink disease virus causes interstitial pneumonia 
in young mink and chronic immunological disorders in 
adult mink, but infection can also occur without clinical 
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signs. Gray fox amdovirus was identified in foxes with ab- 
normal gait and muscle inflammation [38] . In fecal material 
of red fox sample 40, 1 1 sequences were detected that had 
the closest similarity to viruses belonging to the genus 
Amdovirus. Based on these sequences, specific primers 
were designed and sequences of two partial ORFs (left 
ORF, putative non-structural gene 1 and 5' end; 633 nt 
and right ORF, putative viral protein 2; 830 nt, Genbank 
accessions KF823809 and KF823808) of this novel virus, 
tentatively called Red fox fecal amdovirus, were confirmed 
by Sanger sequencing (Figure 5A). Phylogenetic analysis 
and calculation of the pairwise identities of the deduced 
amino acid sequence of the partial VP2 gene revealed that 
this virus is most closely related to Grey fox amdovirus, 
with pairwise identities of respectively 83% on the amino 
acid level and 80% on the nucleotide level (Figure 5B). 



Red fox fecal kobuvirus 

Kobuviruses are small single stranded positive sense RNA 
viruses belonging to the family of Picornaviridae. At 
present, the genus Kobuvirus consists of three species, 
Aichivirus A, B and C. Aichivirus A, B and C were de- 
tected in enteric samples from humans, cows and pigs re- 
spectively [39-41]. In addition, kobuviruses were detected 
in fecal samples of a number of other species recendy 
[42]. Although seroprevalence studies suggest that infec- 
tions with kobuviruses are common, their role as a pri- 
mary pathogen is unknown [43,44] . 

Sequences with the closest similarity to kobuviruses 
were detected in a common genet (sample 13) and in 
two red foxes (samples 37 and 39) (Table 1). Kobuvirus 
sequences detected in the common genet had the closest 
similarity to kobuviruses detected in mice and sewage 
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Figure 5 Phylogenetic analysis of Red fox fecal amdovirus. 

Overview of obtained viral sequences of Red fox fecal amdovirus 
using Aleutian mink disease virus as a reference (A). Phylogenetic 
neighbor-joining tree with p-distance and 1,000 bootstrap 
replicates of the deduced amino acid sequence of the partial VP2 
gene of Red fox fecal amdovirus and VP2 genes of various other 
viruses of the genus Amdovirus detected in mink and gray foxes 
and viruses of the genus Parvovirus detected in red and blue 
foxes and dogs (B). Indicated are Genbank accession numbers 
and bootstrap values. 



[7,45], while kobuvirus sequences detected in the foxes 
were nearly identical to kobuviruses detected in healthy 
and diarrheic dogs [43,46,47]. Using primers described 
previously [46], the partial 3D region of a variant kobu- 
virus detected in feces of a red fox (sample 37), tenta- 
tively called Red fox fecal kobuvirus (Genbank accession 
KF823813) was amplified and sequenced. Alignment and 
phylogenetic analysis of this Red fox fecal kobuvirus with 
various other kobuviruses indeed confirmed that detected 
sequences were very similar to canine kobuviruses, with 
an identity of 97% on the nucleotide and 100% on the de- 
duced amino acid level with Canine kobuviruses detected 
in Italy recently [46] (Figure 6). The high similarity be- 
tween kobuviruses detected in dogs and foxes suggest that 
kobuviruses of dogs and red foxes can be easily transmit- 
ted from one host to another or that transmission of these 
viruses has occurred relatively recently. 



Novel picobirnaviruses 

Picobirnaviruses are small, non-enveloped, bisegmented 
double-stranded RNA viruses. These viruses have been 
often detected in fecal samples of humans and various 
animal species with and without disease [7,14,48,49]. In 
the present study, sequences that had the closest similar- 
ity to viruses belonging to the family Picobirnaviridae 
were detected in fecal samples collected from two com- 
mon genets, two European minks and two red foxes 
(Table 1). Based on 454-sequencing reads, the (partial) 
coding sequences of the RNA-dependent RNA polymer- 
ase (RdRp) gene of a novel picobirnavirus detected in 
one common genet (sample 14; 811 nt, Genet fecal pico- 
birnavirus, Genbank accession KF823812) and of two 
novel picobirnaviruses detected in a red fox (sample 40, 
Red fox fecal picobirnavirus 40-1; 1560 nt and 40-2; 
1669 nt, Genbank accessions KF823810 and KF823811) 
were obtained (Figure 7A). Alignment and phylogenetic 
analysis of these viruses with RdRp gene sequences of 
similar length of other picobirnaviruses showed that ob- 
tained sequences of the RdRp gene were most closely re- 
lated to picobirnaviruses detected previously in rodents 
and otarines (Figure 7B). 

Various sequences of known and novel viruses were 
identified but no known zoonotic viruses were detected. 
In addition to the viral sequences characterized in the 
present manuscript, a novel rhabdovirus was discovered 
in feces of a red fox (nr 40, Bodewes et al., manuscript 
submitted for publication), which was divergent from 
currently known rhabdoviruses, with pairwise identities 
on the amino acid level of the L gene of this novel rhab- 
dovirus with other viruses of the family Rhabdoviridae 
of only 35% or less. 

The presence and detection of the obtained viral se- 
quences is potentially biased by multiple factors, includ- 
ing ecological factors, virus concentrations, a PCR bias, 
a bioinformatics analysis bias and possible contamin- 
ation of laboratory kits since we also identified se- 
quences with the closest similarity with NIH-CQV/PHV, 
although only in one sample [16]. Since for almost all 
animals only fecal material was available, it is unknown 
whether these novel viruses have caused disease in the 
host or if they might potentially transmit to domestic 
animals or humans. Additional studies based on collec- 
tion of blood samples and complete necropsies are ne- 
cessary to elucidate the pathogenicity and epidemiology 
of these viruses. Of interest, a number of species evalu- 
ated in the present study have been the focus of previous 
viral metagenomics studies in the Netherlands [9,14]. 
Still, high numbers of previously unknown viruses were 
detected. This might be due to the relatively low number 
of animals sampled in each study or due to the geo- 
graphic distance between the Netherlands and Spain. 
Further studies using samples of more countries are 
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necessary to elucidate the geographical spread of each of 
these viruses among animals in Europe and beyond. In 
addition, in spite of their relevance as potential carriers 
of pathogens, only a few recent studies have conducted 
thorough metagenomic analysis of fecal material on wild 
carnivore species [9,14,50]. 

Conclusions 

The present study highlights the viral diversity that is 
present in fecal material of wild carnivores. Results ob- 
tained in this study are a step forward in order to improve 
the limited knowledge about virus diversity present in wild 
carnivores in Europe and will help to get a better under- 
standing of potential cross-species transmission of viruses 
between diverse hosts [51] and provide information about 
potential pathogens [52]. 

Materials and methods 

Sample collection 

Rectal swabs or fecal specimens were collected from 42 
carnivores in the Basque Country and La Rioja regions, 



Spain (Figure 1; Table 1). Samples from the Canidae, 
Mustelidae, Viverridae and Felidae families of the order 
Carnivora, including American mink (n = 10), Common 
genet (n = 7), Eurasian badger (n = 4), Eurasian otter (n = 1), 
European mink (n = 9), European pine marten (n = 2); 
European polecat (n = 3), red fox (n = 4), stone marten 
(n = 1), and European wild cat (n = 1) were analysed in 
the present study. Carnivore samples were obtained from 
found-dead animals, (road casualties and/or poaching) col- 
lected by authorized veterinarian personnel of the Wildlife 
Rehabilitation Centre of Martioda (Alava Regional Council. 
Department of Environment; Biodiversity section). Add- 
itionally, samples from live trapped individuals were ob- 
tained in the framework of a systematic population study 
of the Endangered European mink {Mustela lutreola) and 
the control and eradication program of the invasive alien 
American mink in Spain, for other purposes than this 
study (Ref: 3088046-CONSERVISON. Conservation actions 
for the Endangered European mink: National strategy de- 
velopment, ex-situ conservation program and American 
mink eradication. TRAGSATEC - Spanish Ministry of 
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Agriculture, Food, and Environment). American minks 
were humanely killed according to ethical and legal re- 
quirements and in agreement with the Spanish law in 
order to control alien invasive species (Law 42/2007 and 
Real Decree 1628/2011). All the trapping, handling, 
sampling and American mink culling was performed by 
specialized wildlife researchers and veterinarians with 
the permission of regional wildlife authorities (La Rioja 
Government. Agriculture, Livestock and Environmental 
Council. General Direction of Natural Environment. 
Nature Conservation and Planning Service: Ref: LL/aic; 
Alava Regional Council. Department of Environment. 
Biodiversity Section: Ref: 13/32) and in line with the 
laws and ethical protocols governing wildlife manage- 
ment (Law 42/2007; Decree 14/2002; Foral Order 180/2003 
and 322/2003). No animals were sacrificed for the pur- 
poses of this study. Therefore, a formal approval by an In- 
stitutional Animal Care and Use Committee was not 
necessary. Following collection, samples were directly 
stored at -20°C and were stored at -70°C within 2 months 
after collection until further processing. 



Sequence independent RNA and DNA virus screening of 
collected samples 

Samples were processed for viral metagenomics as de- 
scribed previously [9,49]. In brief, samples were depleted 
from host nucleic acids and filtered through a 0.45 uM fil- 
ter. Subsequently, RNA and DNA were extracted using 
the Nucleospin RNA XS kit (Macherey-Nagel) and the 
High Pure viral nucleic acids kit (Roche). First and second 
strand synthesis and random PCR amplification were per- 
formed. PCR products were purified and processed for 
next-generation sequencing with a 454 GS Junior Instru- 
ment (Roche). Obtained reads were assembled using de 
novo assembly in CLC Genomics Workbench 5 (CLC Bio) 
and contigs and individual reads were analyzed by BLASTN 
and BLASTX respectively. Cut off E-values for significant 
virus hits for BLASTN and BLASTX were respectively 
1.0 x 10~ and 1.0 x 10~ . Based on the taxonomic origin 
of the best-hit sequence, classification of the sequences 
was performed in MEGAN 4.70.4 [53]. Obtained reads 
were deposited at the European Nucleotide Archive under 
archive number PRJEB4910. 
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PCR amplification and sequencing 

Based on obtained next-generation sequencing data, spe- 
cific primers of the novel rhabdovirus, theilovirus, and 
amdovirus were designed to confirm and partially extend 
the obtained sequence data. Products were sequenced as 
described previously and primer sequences are available 
upon request [49]. 

Phylogenetic and protein analysis 

Nucleotide and/or deduced amino acid sequences of the 
novel kobuvirus, theilovirus, amdovirus, and picobirna- 
viruses were aligned using ClustalW in MEGA5 with de- 
fault parameters [54], while deduced amino acid sequences 
of the novel phleboviruses were aligned using MAFFT (vs 
7; http://mafft.cbrc.jp/alignment/software/) with the E-INS- 
I algorithm and otherwise default parameters. Phylogen- 
etic analysis of the novel theilovirus, kobuvirus, amdovirus 
and picobirnaviruses were performed using a neighbor- 
joining tree with 1000 bootstrap replicates in MEGA5, 
while phylogenetic analysis of GP and NP genes of the 
novel phleboviruses was performed using a maximum- 
likelihood tree (WAG + F + I + G model) with 100 boot- 
strap replicates. Alignment of the amino acids of the 
major surface structures of the theiloviruses was per- 
formed with JalView version 2 [55]. Prediction of cleav- 
age sites of picornaviral proteases was performed using 
NetPicoRNA 1.0 (www.cbs.dtu.dk/services/NetPicoRNA/). 
Other conserved motifs were predicted with MEME 4.9.0 
(http://meme.nbcr.net/meme/). The presence of conserved 
domains was evaluated using the Conserved Domains 
Database [56]. 
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